Strong Correlations Produce the Curie- Weiss Phase of Na^CoC^ 
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Within the t-J model we study several experimentally accessible properties of the 2D-triangular 
lattice system Na^CoC^, using a numerically exact canonical ensemble study of 12 to 18 site tri- 
angular toroidal clusters as well as the icosahedron. Focusing on the doping regime of x ~ 0.7, we 
study the temperature dependent specific heat, magnetic susceptibility and the dynamic Hall coef- 
ficient Rh (T, uj) as well as the magnetic field dependent thermopower. We find a crossover between 
two phases near x ~ 0.75 in susceptibility and field suppression of the thermopower arising from 
strong correlations. An interesting connection is found between the temperature dependence of the 
diamagnetic susceptibility and the Hall-coefficient. We predict a large thermopower enhancement, 
arising from transport corrections to the Heikes-Mott formula, in a model situation where the sign 
of hopping is reversed from that applicable to Na^CoC^. 

PACS numbers: 72.15.Jf, 65.90.-K, 71.27.+a 
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The physics of the 2-dimensional triangular lattice sys- 
tem sodium cobaltate Na^CoC^ (NCO) is fascinating[l|, 
combining strong electron correlations and thermoelec- 
tric physics. A Curie- Weiss metallic phase for dopings 
x ~ 0.7 has been reported @ where the observed physical 
variables display an unusual mix of behaviors that are 
hybrid between those of good metals and of insulating 
systems challenging theory severely. The thermopower 
of NCO is nearly ten times higher than expected from 
typical metals; generating excitement in the engineering 
and material science communities, particularly applied to 
thermoelectric devices. Here we show that strong elec- 
tron correlations, along with the geometrically frustrated 
lattice of NCO, hold the key to explaining this mysterious 
state of matter. We examine several experimentally ac- 
cessible properties of NCO within the t-J modelU0,0. 
We find that strong electron correlations capture the es- 
sential physics and our results compare well with exper- 
iment. 

In NCO the low spin Cobalt ion valence fluctuates be- 
tween a Co 4+ (spin 1/2) and a Co 3+ (spin 0) configura- 
tion; the number of Co 3+ states is precisely x. The Co 
ions form a triangular lattice, and photoemission0, Q 
is consistent with a single, hole-like, band with hopping 
t < and n = 1 + x electrons satisfying the Luttinger 
volume count. 

The t-J model describes strongly correlated electron 
systems by forbidding double occupancy of lattice sites. 
We apply this model to NCO after an electron-hole trans- 
formation, requiring t — > —t and hole doping |1 — n\. A 
non-zero J couples nearest neighbor electrons via their 
spin degree of freedom. For such strongly correlated sys- 
tems perturbation theory is doomed to failure from the 
outset and we make progress through numerical exact di- 
agonalization on systems containing 12, 14, and 18 sites 
on toroidal clusters (periodic boundary conditions (BC)) 
and on ladder clusters (open BC in one direction) Q • 
Thermodynamics is considered within the canonical en- 



semble. 

The Hilbert spaces of these finite systems are very big 
(up to ~ 80,000 states) and grow exponentially with the 
number of sites. Therefore, all available symmetries are 
used to reduce the dimension of the matrices that arise 
to large but manageable proportions. However, Peierls 
phase factors are needed to describe an applied mag- 
netic field, which breaks or reduces the translational in- 
variance, thereby limiting us somewhat. By using a ju- 
dicious choice of the BC and of phases on bonds [lOJ, we 
achieve a fairly small non-zero flux per plaquette of n/Nf, 
where Nf is the total number of triangular faces on the 
lattice^. The ladder systems, however, enable an in- 
finitesimal flux to be chosen. 

For NCO, photoemission supports a value for the hop- 
ping of t = —100 K and we adopt it in this work. This 
value is suggested by the ARPES data 6] on the loss of 
coherence of the quasiparticles as well as the dispersion 
in the composition range x ~ .7. The T dependence of 
the chemical potential n(T) — /i(0) is another route to 
estimating t [llj. 

Fig. 1(a) shows the electronic specific heat C V (T), and 
is compared with that for non-interacting electrons with 
the same hopping. We find that the effect of correla- 
tions is a shift in the peak to a smaller temperature and 
suppression of its overall weight. This is expected since 
the Gutzwiller projection in the t-J model reduces the 
number of available states and hence the entropy. Due 
to a finite system size induced gap in the spectrum, we 
expect an exponential behavior of C V (T) for T < 20 K, 
the typical gap value. Taking this into account, we are 
able to extract the linear electronic contribution jT. The 
value of 7 is enhanced by ~ 1.5 over the non- interacting 
value; this enhancement depends only slightly on J (ne- 
glecting the exponential increase at T < 20K due to the 
finite-system-induced gap) and varies with system size by 
0.2. J = 40 K (i.e., 0.4|i|), which is fixed by the experi- 
mental system through a comparison of the Curie- Weiss 
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FIG. 1: (color online) Specific heat and susceptibility, (a), 
specific heat C V (T) for x = 0.72, computed on the 18-site 
cluster, comparison of J = (bottom) with J/\t\ = 0.4 (mid- 
dle) and bare (Hubbard U = 0) specific heat (top), dotted 
straight lines show linear fits and 7 values for J = and 
U — in units of mJ/(mol-K 2 ). (b)-(f), susceptibility x(T) 
for dopings around x ~ 0.7. The dotted curves indicate the 
bare susceptibility, and arrows indicate the evolution of J/\t\ 
from to 0.5 in steps of 0.1 (red to yellow). Note the change 
of scale in different panels. These results combine two dif- 
ferent clusters, a 12-site torus [x = 0.58, 0.67, 0.75) and a 
14-site torus (as = 0.71, 0.79). The difference in x = 0.71 and 
x — 0.75 shows that x(T) transitions from a Curie- Weiss to 
Pauli paramagnetic behavior in this range. 



temperature with computations0,|i|]. 

In Fig. l(b)-(f) we present the spin susceptibility 
x(T) for several dopings around x ~ 0.7. In the band 
limit x — > 1, as in the upper two panels (as well as 
results not shown), we find the expected weakly T de- 
pendent but J insensitive Pauli paramagnetic behavior. 
When x is lowered below x — 0.75 (bottom three panels), 
x(T) shows strong Curie- Weiss-like T and J dependence, 
and is significantly renormalised from the non-interacting 
value at low T. This indicates a crossover to the strong- 
correlation induced local moment behavior for x < 0.75 
which closely parallels experimental findings 0. In this 

Curie- Weiss phase, the behavior at high T is described 

2 2 

by the Curie- Weiss form x(T) = | ^ w ^ n a ne g a_ 

tive Weiss temperature 9 and effective magnetic moment 
Peff, v — V/N is the unit cell volume. When continuing 
the analysis to 1 — > 0, antiferromagnetic (AFM) corre- 
lations increase and we find that 9{x,J) — —cJ e ft(x) 



where J e ff(x) = J(l + c"x\t\) + c'x\t\, with c = 4.0, 
d = 0.01425, and c" = -0.9175. The d term origi- 
nates in the kinetic antiferromagnetism of the frustrated 
lattice^]], an d signifies that even in the absence of J, 
there is a tendency for AFM order, i.e., in a direction 
opposite to the usual Nagaoka mechanism for the square 
lattice [l^. 

Experimentally, the Hall coefficient of NCO is re- 
markable in many respects. Most striking is the un- 
bounded linear increase with temperature of the Hall co- 
efficient Rh- To understand this we perform the brute 
force exact summations of Kubo's formulae for various 
conductivities by introducing a level width, i.e., a 
broadening lu — > lu + ir\ with 77 equal to the mean en- 
ergy level spacing. In addition, we evaluate the high fre- 
quency limit 3] of Rh (called R* H ) for all T. Recall that 
the high T estimates of R* H led to a prediction^, 0] of 
the linear T dependence of the Hall constant for NCO, 
which was successfully verified^^. We are thus able to 
provide a purely theoretical benchmarking of this idea as 
well, subject of course to the limitations of the finite size 
clusters. 

Focusing on the region of doping around x ~ 0.7 Fig. 
2(a) shows the Hall coefficient as a function of temper- 
ature and frequency. We find that the Hall coefficient 
is relatively insensitive to freq uency, in keeping with the 
original expectations 0, 0, Il7| . All curves show a mini- 
mum near T = 100 K, and an unbounded linear increase 
for T > 200 K. The slope of R* H as found in the clus- 
ters is in agreement with results from high-temperature 
expansions [3j. The experimental curve, unlike theory, 
has a change of sign and also a pronounced minimum 
at T ~ 100 K. We cannot reproduce the change of sign. 
While we could fit the high T slope more accurately, it 
requires a smaller value of hopping t ~ 35 K, as already 
noted Q. Such a choice would make most other vari- 
able fits less sensible. Hence we conclude that while the 
data^ is largely as expected from theory, in detail it is 
still not possible to reconcile the change in sign as well 
as the magnitude of the slope with theory. Therefore, 
the Hall coefficient still offers a considerable challenge to 
both theorists and experimentalists. 

In an alternate effort to understand further the 
data^^ at lower T ~ 100 K, we note that the curvature 
of the (Landau) diamagnetic susceptibility \d (obtained 
by inserting Peierls phases) and the Hall constant are 
curiously related in our computation via 



Td 2 Xd /dT 2 = f(x)d 2 R H /dT 2 



(1) 



where f(x) is a function depending on doping x. We in- 
tegrate equation (JTJ to arrive at the Hall constant from 
our Xdi yielding a good overall fit of data. The problem 
with the slope and the negative intercept at low T are for- 
given in this approach since one fits the two constants, 
but it does capture the ubiquitous minimum at T ~ 100 
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FIG. 2: (color online) Hall coefficient and thermopower. (a), 
Comparison of several results for the Hall coefficient Rh(T) 
at x = 0.75 with experiment in ref. |lql at x = 0.71 (red 
squares): R* H (blue dot-dash), Rh(u) = 0) (blue dashed), Rh 
(orange solid) derived from \d (at x = 0.83, ladder), the 
dc-limit required a broadening of the frequency u — > u) + if] 
with 77 « 3|t| to eliminate finite-size artifacts. All results 
are for 12-site clusters and J/\t\ = 0. (b), Infinite frequency 
thermopower S* vs. T for a 12-site torus at x = 0.75 and 
x = 0.67. The solid black and dashed blue lines correspond 
to J/\t\ = 0, and 0.4 at x = 0.75, respectively, while the solid 
orange and dashed-dotted red lines correspond to J/\t\ = 0, 
and 0.4 at as = 0.67, respectively. S*(T) for t = -100 K rel- 
evant for NCO. The diamonds and stars represent measured 
thermopower for NCO at x = 0.68 from ref. and ref. [l3j. 
(c), Our prediction for S*(T) for the case when the sign of 
the hopping is reversed (t = 100 K.) 



K, with a slight J dependence |3- While the Landau dia- 
magnetism is very small compared to the paramagnetic 
susceptibility for narrow band systems, it might yet be 
accessible to experiments. This is so since it is anisotropic 
in space (being planar), in contrast to the isotropic Pauli 
susceptibility, and hence torque magnetometry 18] could 
help disentangle these terms. 

The thermopower of NCO is striking in its large 
magnitude [l^| ~ 100^y/K and also in its surprising sen- 
sitivity to an applied magnetic fieldQ. A formulation 
for the thermopower has been recently given |l7^ in the 
high frequency limit in the same spirit as the Hall con- 
stant. We note that the thermopower can be written 
as a sum of a transport and the Heikes-Mott term as 
S(oj,T) = S Tr (u,T) + Shm(T), where the Heikes-Mott 



term is Shm(T) — q y , and the transport term is 

^■toSM^. M^-fW . (2) 

Here q e — — |e| is the electronic c harg e and f xx is the 
diagonal part of the stress tensor [jl 1 17j . The term Shm 
is entropic in original. Detailed expressions for A as an 
expectation of a many-body operator in the t-J model 
are given in equation (83) of ref. (l7|. and in a longer 
workQ- 

The main approximation in ref. 17] is to use the high 
u) limit of S(ui,T) (called S*), and is expected to be nu- 
merically quite reasonable, in parallel to the behavior of 
the Hall constant reported in this work. At low T the two 
contributions to S vanish separately as we have written 
them, and in general we find that for t < (the case 
of NCO) the entropic part is by far the dominant term. 
For the opposite case (t > 0) the transport term comes 
into play in a dominant way, and leads to very interesting 
behavior as we shall show. 

We first discuss the thermopower in the absence of 
magnetic field. Fig. 2(b) shows the T and J depen- 
dence of S* for x relevant to the Curie- Weiss phase. The 
Heikes-Mott term dominates over the transport term and 
the frequency dependence of the thermopower (evaluated 
via Kubo formulae) is found to be quite weak in this 
range of doping for t < relevant to NCO 8], and thus 
the high frequency approximation is as good as exact. 
We find our results compare well with the experimen- 
tal data0, El- We now sharpen the prediction (from 
Eq. (88) of ref. 0) that reversing the sign of hopping 
leads to a maximum in S* as a function of T. Fig. 2(c) 
indeed shows such a behavior and provides an estimate 
of the expected enhancement in S* , a factor of 2 to 3 
over the frustrated case of NCO and also of the expected 
temperature scale T max ~ \t\. It would be interesting to 
check this experimentally. A part of the enhancement 
for t > arises from the prominent peak in the single 
particle density of states coming close to the fermi level. 
Our numerics show that interactions further amplify this 
effect considerably (a factor of ~ 2). 

Next we consider the effects of a magnetic field of 
strength B (perpendicular to the plane) on the ther- 
mopower. The field induced change of the transport term 
arises from the Peierls factors, and is found to be a very 
small fraction (~ 0.001) of the change in the Heikes-Mott 
term. The field dependence of the chemical potential is 
significant and responsible for the overall change. In Fig. 
3(a) we show the normalised Shm(B, T) for several dop- 
ings as function of B/T. This scaling with B/T is very 
similar to the one found in ref. |2| , both qualitatively and 
quantitatively. In Fig. 3(b)-(d) we find a crossover for 
x ~ 0.75 from a weakly independent S(B, T) at x > 0.75 
to the Curie- Weiss phase where S is greatly suppressed 
by B. This crossover is similar to the one seen in the 
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FIG. 3: (color online) Magnetic field dependence of the 
thermopower. (a), Shm(B) normalised to Shm{0) shows 
the relative suppression of Shm with B for temperatures 
T = (100, 120, 140, 160, 180,200) K (blue to red), horizontal 
axis is rescaled by 1/T as done in ref. 0. (b)-(d), Shm(T) 
for several doping values around x — 0.71 and several values 
of B from T to 13 T (red to orange). Note that for x > 0.75 
the field-dependence becomes very weak. The insets display 
the full temperature behavior of Shm out to T = 400 K. All 
data for J = 0. Finite J reduces the field suppression. 



T dependence of the spin-susceptibility. These results 
confirm the interpretation 2] in terms of spin-entropy as 
the leading contribution to the field-suppression and also 
provide a guide to what one can expect at high magnetic 
fields that are not accessible experimentally. 

We next make a few comments about the overall be- 
havior of the thermopower and Hall constant as a func- 
tion of filling in strong coupling theories. As noted 
earlier 0, H(J, while band theory predicts a single zero 
crossing of the Hall constant as n increases from 0^2, 
a strongly correlated electron system has three cross- 
ings. In precisely the same sense the thermopower also 
has three zero crossings for a correlated band system. 
We demonstrate this for the triangular lattice. In Fig. 
4(a)-(b) we show the filling dependence of both the Hall 
constant and the thermopower. The divergence of the 
Hall constant at half filling is forced by the carrier freeze 
out accompanying the Mott-Hubbard gap, and is less 
pronounced in the triangular lattice than in the square 
lattice|9|. For the thermopower we expect similar results 
for the square lattice. 

The hopping parameter t (~ 100 K) used here to recon- 
cile a considerable body of experimental data with theory, 
is much smaller than that in high T c systems t ~ 3000 
K. This value is also smaller, by a factor of 10, than the 
one found in LDA studies |2lL |22j , and should motivate 
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FIG. 4: (color online) Hall coefficient and thermopower as a 
function of filling n. (a), Filling dependence of R*h(T = 0) for 
a 12-site torus (black circles) and icosahedron (red diamonds) 
as well as Rh(T — 0,lj = 0) for the 12-site torus (blue trian- 
gles) all for J/\t\ = 0. (b), Filling dependence of S*(T) for 
T = 100 K, 200 K, and 1000 K (blue triangle, red star, black 
triangle) for J/\t\ — on a 12-site torus. The solid black line 
is Shm at infinite T. The dotted lines in both panels are 
guides to the eye indicating zero Rh and S and half filling 
(n = l). 

further studies of the difficult problem of projecting a 
multi-band system onto a single band model systemati- 
cally. In conclusion, our work shows that strong correla- 
tions account for the observed dramatic behavior of the 
Curie- Weiss metallic phase in considerable quantitative 
detail. 

We thank N. P. Ong and Y. Wang for valuable 
discussions. This work is supported by Grant No. 
nsf-dmr0408247. 
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